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ABSTRACT 

We have investigated optical properties of 4248 infrared galaxies (IRGs) by positionally 
matching data from the Infrared Astronomical Satellite (IRAS) and the Sloan Digital Sky Sur- 
vey Data Release 3 (SDSS DR3). Due to the large sky area coverage of both surveys, we have 
obtained an opportunity to study an unprecedentedly large number of 4248 infrared galax 
ies spanning four orders of infrared luminosity range 1O 9 L sC Lg_ iooo>m ^ 
Our sample includes three hyper luminous infrared galaxies (ig-ioooum 
large number of 178 ultra luminous infrared galaxies (10 12 L Q < L 



1O 1357 L . 
> lO i3 L ) and the 



8-1000>m 



1 INTRODUCTION 



One of the most important achievements of the infrared astronom- 
ical satellite (IRAS; Neugebauer et al. 1984) is the discovery of 
a significant population of luminous infrared galaxies (Soifer et al. 
1984; Sanders & Mirabel 1996; Sanders 2003). Especially, the most 
energetic objects, the ultra luminous infrared galaxies, has been 
a subject of extensive debate (ULIRGs; Ls-iooo/im > W 12 Lq; 
e.g., Sanders et al. 1988; Borne et al. 2000; Kim et al. 2002) since 
ULIRG's bolometric luminosity is comparable to that of QSOs 
(Soifer et al. 1987) and dominate the top end of the galaxy luminos- 
ity function. The high infrared (IR) luminosity is attributed to dust 
emission, with the heating source being varying combinations of 
starbursts (Genzel et al. 1998; Lutz et al. 1998) and active galactic 
nuclei (AGNs; Veilleux et al. 1999a,b;Nagar et al. 2003). The lumi- 
nosity functions of infrared galaxies show that at Li r > 10 Lq, 
infrared selected galaxies become more numerous than optically 
selected starburst and Seyfert galaxies of comparable bolometric 
luminosity (Soifer et al. 1987). At L ir > 10 12 L©, ULIRGs ex- 
ceed the space densities of QSOs by a factor of 1.5-2 (Sanders et 
al. 1989). Some investigators have proposed that ULIRGs may rep- 



< 10 1J L Q ). 

In addition, it is important to have statistical number of lower luminosity infrared galaxies 
(is-iooo/jm ~ 10 10 I/q) in order to link the well-studied ultra luminous infrared galaxies to 
normal star-forming galaxies. Our findings are as follows: (i) we found that more IR luminous 
galaxies tend to have smaller local galaxy density, being consistent with the picture where lu- 
minous IRGs are created by the merger-interaction of galaxies that happens more often in 
lower density regions; (ii) the fractions of AGNs increase as a function of Li r ; (iii) there is a 
good correlation between L.- Lr and Lipiii] for AGNs, suggesting both of the parameters can 
be a good estimator of the total power of AGNs; (iv) a good correlation is found between 
Li r and optically estimated star formation rate (SFR) for star-forming galaxies, suggesting 
Li r is a good indicator of galaxy SFR. However, caution is needed when SFR is estimated 
using L^, i-e., high SFR galaxies selected by Li r is frequently merger/interaction, whereas 
high SFR galaxies selected by optical emission is often normal spiral or Magellanic-cloud 
like irregular galaxies, (v) more IR luminous galaxies have slightly larger Ha/H(3 ratio; (vi) 
more IR luminous galaxies have more centrally-concentrated morphology, being consistent 
with the morphological appearance of galaxy-galaxy merger remnants. Optical image of ul- 
tra/very luminous infrared galaxies also show frequent signs of merger/interaction, (vii) com- 
parison with the SED synthesis models indicate that majority of luminous infrared galaxies 
(-Lg_iooo>m > 1O 11 L ) may be in a post-starburst phase, sharing a similar (but not the same) 
merger/interaction origin with E+A galaxies. 

resent the initial, dust-enshrouded stages of QSOs (e.g., Sanders et 
al.1999). Some ULIRGs (especially cool ULIRGs) also appear to 
be forming moderately massive (L*) field elliptical galaxies (e.g., 
Melnick & Mirabel 1990;Kormendy & Sanders 1992; Genzel et 
al. 2001). Their near-IR light distributions often fits an r^-law 
(Wright et al. 1990; Stanford & Bushouse 1991; Doyon et al. 1994; 
Scoville et al. 2000; Veilleux et al. 2002). Based on both optical 
and infrared morphology, many studies found that the majority of 
ULIRGs are strongly interacting or advanced merger systems (e.g. 
Borne et al. 2000). 



On the other hand, lower luminosity infrared galaxies 
(lO 9 I/0 < Li r < W 12 Lq) have not been the subject of much 
scrutiny (however, see Kim et al. 1995; Wu et al. 1998a,b; Arribas 
et al. 2004). However, lower luminosity infrared galaxies are also 
an interesting class of galaxies for various reasons: it has been re- 
vealed that the properties of ULIRGs appear to correlate with IR 
luminosity (Sanders & Mirabel 1996). Therefore, in order to under- 
stand infrared galaxies (IRGs) as a whole, we need to study IRGs 
at a lower luminosity range. The luminosity of lower luminosity 
IRGs investigated here are between ULIRGs and more popular spi- 
ral galaxies (W 9 Lq < Li r < 10 12 Lq ). Therefore, lower luminos- 
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ity IRGs might provide a link between normal spiral galaxies and 
infrared luminous galaxies. In addition, lower luminosity IRGs are 
much more frequent than ULIRGs, and thus, it is important to study 
VLIRGs in order to understand cosmological evolution of galaxies 
as a whole. 

However, previously, it has been difficult to study statistical 
number of lower luminosity IRGs. It has been simply difficult to 
obtain telescope time to follow-up thousands of IRGs both spec- 
troscopically and photometrically. One of the pioneers is Arribas 
et al. (2004) who observed 19 very luminous infrared systems 
(VLIRGs; 10 n L© < L ir < 10 12 L o ) in order to probe an im- 
portant link between ULIRGs and VLIRGs. They found that the 
photometric properties of ULIRG and VLIRG samples, considered 
as a whole, was indistinguishable at optical wavelengths. Wu et al. 
(1998a,b) carried out spectroscopic observation of 73 VLIRGs to 
find that 60% of VLIRGs exhibit AGN-like spectra, and that 56% 
of VLIRGs show strong merging/interaction signatures. The frac- 
tions are smaller, but these are similar characteristics to ULIRGs. 
With the advent of the Sloan Digital Sky Survey (SDSS), we have 
a chance to step forward. By matching up ~360,000 spectroscopic 
galaxies in the SDSS Data Release 3 (SDSS DR3; Abajazian et 
al. 2004) with the IRAS data, we have obtained the opportunity to 
study optical photometric and spectroscopic properties of unprece- 
dented number of 4248 infrared galaxies, including three hyper lu- 
minous infrared galaxies and 178 ultra luminous infrared galaxies. 
Unless otherwise stated, we adopt the best-fit WMAP cosmology: 
(h, n m , fix,) = (0.71, 0.27, 0.73) (Bennett et al. 2003). 



2 DATA 

2.1 Matching up the IRAS and the SDSS catalogs 

In order to create a large sample of galaxies with both optically 
measured redshift and mid- and far-infrared photometry, we have 
positionally matched up the IRAS galaxies (IRAS faint source cat- 
alog Version 2; Moshir et al. 1992; FSC92 hereafter) and the SDSS 
galaxies (The Sloan Digital Sky Survey Data Release 3; Abazajian 
et al. 2004) with spectroscopically measured redshifts. Moshir et 
al. (1992) provided more accurate fluxes and positions than that 
of 1988 and is considered as the standard reference catalog for 
IRAS sources (However, note possible inaccuracies in using cat- 
alog line fluxes). From IRAS catalog, we exclude IRAS sources 
with /so < /12 to remove contamination from infrared bright 
stars which are known to peak around 12 iim (Cohen et al. 1987). 
This criterion has been known to be efficient in selecting galax- 
ies, while excluding most Galactic objects. ULIRGs typically have 
/60//12 > 20. Positional accuracy of the IRAS catalog in cross- 
scan direction is ~10 arcsec. Therefore, if an IRAS galaxy has a 
counterpart in the SDSS within 20 arcsec, we regard it a match. 
Total number of galaxies with both the IRAS and the SDSS detec- 
tions is 4248, which is one of the largest data-sets with both optical 
redshift and the IR photometry. When we construct a volume lim- 
ited sample, we select a redshift and absolute magnitude range of 
z ^ 0.06 and M r ^ —19.42 to maximize the number of galax- 
ies in the volume limited sample, resulting in 1666 galaxies in the 
sample. 



2.2 Infrared Luminosity 

For each galaxy in the matched catalog, we compute infrared flux 
using the following equation. 
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Figure 1. The histogram of total infrared luminosity, (Li r ). 

F lr = 1.8 x (13.48 x /12+5.I6 x f 25 + 2.58 x / 60 + 1.0 x /u,o)(l) 

The IRAS fluxes are taken directly from the FSC92. Since many 
of the high luminosity galaxies emit a significant portion of their 
infrared luminosity shortward of 40 ttm, this expression provides 
a significantly better determination of the total infrared flux than 
the more commonly used Fi r determined by fitting a single tem- 
perature dust model to the 60 fim and 100 /im fluxes (Kim & 
Sanders 1998). This flux, Fi r , is then converted to luminosity 
(Li r ) using the redshift measured by the SDSS. To keep consis- 
tency with previous work, we use h — 0.75 and qo — only 
for this computation. We show the histogram of total infrared lu- 
minosity (Li r ) in Figure Q We have identified three hyper lumi- 
nous infrared galaxies (HLIRG; 10 13 ' L© < Li r ) and show their 
images and spectra in Figure |2| The most luminous object in the 
sample is SDSSJ091345.5+405628 with L ir = 10 13 ' 57 L Q . Our 
sample also includes 178 ULIRGs, which is one of the largest of 
the kind. We show 9 lowest redshift ULIRGs in Figure [5] The 
corresponding optical spectra are shown in Figure [4] In Figure [5] 
Li r is plotted against redshift. In this work, we call galaxies with 
10 13 '°Z/q ^ £ jr as hyper luminous infrared galaxies (HLIRGs), 



galaxies with 1O 12O L 



< Lir ^ 10 Lq as ultra luminous 



1O 12 L as very luminous infrared galaxies (VLIRGs), galaxies 
with 1O 1O ' 5 L < Lir ^ 10 110 Z/Q as moderately luminous in- 
frared galaxies (MLIRGs),and 10 10 '°L Q < L ir ^ 1O 1O5 L as 
lower luminosity infrared galaxies (LLIRGs), respectively. When 
we mention infrared galaxies (IRGs), the sample includes all the 
classes of infrared galaxies plus galaxies with Lt r ^ 10 1 



3 RESULTS 

3.1 The environment of Infrared Luminous Galaxies 

In this subsection, we investigate the environment of luminous in- 
frared galaxies. It has been often suggested that the ULIRGs are 
caused by the major merger of multiple galaxies (Borne et al. 2000; 




Figure 2. Three HLIRGs (10 13 Z/q < Li r ) identified in this work. The upper panel shows g, r, i-composite images. The corresponding spectra with name 
and redshift are presented in the lower panel. Their luminosity is log(Li r ) =13.05,13.57 and 13.19L©, respectively. 
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Figure 5. Total infrared luminosity measured with the IRAS, (Li T ), is plot- 
ted against redshift measured with the SDSS. 

Taniguchi & Shioya 1998). However, these results are based on the 
eye-based scrutiny of the image/morphology of the central IRGs, 
often with disturbed/dynamical tail signatures. As an independent 
way of measuring the environment of IRGs, we study local galaxy 
density distribution of IRGs. 

We use the volume limited sample (z ^ 0.06 and M r ^ 
— 19.42) in this subsection to ensure neighbors are measured in 



the same absolute magnitude range for each and every galaxy in 
the sample. First, we measure local galaxy density in a fixed phys- 
ical scale. In Figure [5] the local galaxy densities are measured in 
0.5 Mpc (top left panel), 1.5 Mpc (top right panel) and 8 Mpc 
(bottom left panel) of radius in angular direction and ±1000 km 
s in the line of sight direction within the volume limited sample 
(z ^ 0.06 and Mr ^ —19.42). Number of galaxies in each colum- 
nar volume is divided by the physical surface area to be converted 
to local galaxy density. The density computed in a fixed physical 
scale allows us to investigate the environment of IRGs at a spe- 
cific scale. In contrast, in the bottom right panel of Figure [6] the 
density is measured using the distance to the 5th nearest galaxy 
in angular direction and ±1000 km s _1 in the line of sight direc- 
tion within the volume limited sample. Since the distance to the 5th 
nearest scales with galaxy density, this density provides us with a 
scale-free estimate of the environment of IRGs, spanning a much 
wider density range of 4 orders. This parameter is often used in 
previous studies (e.g., Goto et al.2003b,c). In Figure|6| the dotted, 
short-dashed, solid, long-dashed lines denote all galaxies in the vol- 
ume limited sample (regardless of IRAS detection), galaxies with 
VLIRGs, MLIRGs and LLIRGs, respectively. ULIRGs are not in- 
cluded in this analysis due to their small number in the volume 
limited sample. 

It is immediately noted that more infrared luminous galaxies 
have smaller local galaxy density at all scales. For example, in case 
of the 5th nearest density (lower right panel), median density is 
1.11, 1.33, and 1.68 Mpc -2 for VLIRGs, MLIRGs and LLIRGs, 
respectively. Compared with the density distribution of all galax- 
ies in the volume limited sample (dotted line; we regard this as 
the field sample), VLIRGs have significantly lower density, and 
MLIRGs and LLIRGs have higher density than the field sample. 
This is a contrasting result compared with that from optical lumi- 
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Figure 3. Examples of g, r, i-composite images ofULIRGs(10 120 L© < L lr < 10 13 L©). The images are sorted from low to high redshift. Only 9 lowest 
redshift galaxies are shown. The corresponding spectra with name and redshift are presented in Fig|4] (The spectrum of the same galaxy can be found in the 
same column/row panel of Fig[4]) 

v3_2) for k-correction. As is expected from the Li r , median M r is 
brighter for the sample with larger Li r . However, the increase in 
r-band flux is only by a factor of ~1.5, compared to the order of 
~ 2 increase in Li r . The M r of MLIRGs and LLIRGs are even 
comparable to that of all galaxies in the volume limited sample. 

3.2 Starburst/AGN classification 

In the following subsections, we investigate optical properties of 
IRGs in detail. However, there are two types of energy sources 
in infrared radiation: starburst and active galactic nuclei (AGN); 
it is critical to separate these two in order to understand physics 
behind infrared radiation. In this work, we separate AGNs from 
starburst galaxies using the optical emission line ratios as de- 
scribed in Kewley et al. (2001). Figure0shows emission line ratios 



nosity. In optical, more luminous galaxies tend to exist in higher 
density regions — often the most luminous ones sit at the center 
of galaxy clusters (e.g., Goto et al. 2002a,b). Perhaps, the total in- 
frared luminosity (Li r ) reflects star formation rate (SFR) of galax- 
ies rather than the total mass in galaxies. Then, this trend can be 
understood as more star-forming galaxies reside in lower density 
region. If the origin of VLIRGs and MLIRGs is merger/interaction 
of multiple galaxies, it is also naturally explained that they pref- 
erentially exist in lower density regions, where relative velocity of 
galaxies are small, and thus, dynamical merging of galaxies hap- 
pens more frequently. It also should be noted that IRGs do not have 
density peaked at the highest density, suggesting majority of IRGs 
are not in galaxy clusters, but in the field region. 

Using the same volume limited samples, we compute me- 
dian M r of each sample in Table Q We use Blanton et al. (2003, 
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Figure 4. Example spectra of 9 lowest redshift ULIRGs. The spectra are sorted from low redshift. Each spectrum is shifted to the restframe wavelength and 
smoothed using a 20 A box. The corresponding images are shown in Fig. [3] (The image of the same galaxy can be found in the same column/row panel of 
Fig|5J. 



Table 1. Median absolute magnitudes in the SDSS r-band (M r ) within the 
volume limited sample (z $C 0.06 and Mr sj —19.42). 



Sample 


Median M r 


N galaxies 


ALL 


-20.10 


20753 


VLIRGs 


-20.61 


289 


MLIRGs 


-20.34 


924 


LLIRGs 


-20.05 


451 



of [onrj(5008A)/H/3 versus [NII](6585A)/Hq. Here, we use the 
emission line flux measured by the SDSS pipeline (Stoughton et al. 
2002). The contour shows distribution of all galaxies in the SDSS 
DR3 with r < 17.77 (regardless of IR detection). The magenta, 
blue, green, and red dots are for ULIRGs, VLIRGs, MLIRGs, and 
LLIRGs, respectively. The dotted line is the criterion used to divide 
AGNs and starbursts proposed by Kewley et al. (2001). We take a 
conservative approach in selecting AGNs. We only regard a galaxy 
as an AGN when all four lines are detected in emission and the line 
ratios are above the dotted line in Fig0 The rest of the galaxies 
(below the dotted lines or with unmeasurable emission lines) is re- 
garded as starburst galaxies. In our sample, there are 810 AGNs and 



3438 star-forming galaxies with the IRAS detections. Both of these 
are one of the largest number of IRGs observed both in optical and 
IR. 

In Figure Q it is also found that more ULIRGs and VLIRGs 
(magenta and blue) are above the dotted line (AGN/starburst cri- 
terion), suggesting more ULIRGs/VLIRGs are AGNs than lower 
luminosity IRGs. To clarify, in Figure|8| we plot fractions of AGNs 
selected in Figure Q as a function of total infrared luminosity. It is 
understood that at a higher luminosity, more galaxies are powered 
by AGNs, especially at Lj r > 10 12 Z/q. In Figure|5| we show frac- 
tions of AGNs as a function of /25//60 IR color. As previously 
suggested (e.g., Veilleux et al. 1999), at /25//60 < 0.2, the frac- 
tions of AGNs decreases, whereas at /25//60 > 0.2 the fractions 
of AGNs are almost constant. 

Previously, it has been suggested that IR galaxies have lower 
excitation than optically selected starburst galaxies (Sarazin 1976; 
Allen et al. 1985; Leech et al. 1989). However, in Figure [7] the 
distributions of dots (IRGs) are consistent with that of the contours 
(all SDSS galaxies). Perhaps, the comparison of excitation ratio 
depends on how you select a starburst sample. Our data show that 
when compared with all the SDSS galaxies, the excitation of IRGs 
is not significantly lower than that of the optically selected galaxies. 
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Figure 6. Normalized distributions of local galaxy density. The local galaxy densities are measured in 0.5 Mpc (top left panel), 1.5 Mpc (top right panel) and 
8 Mpc (bottom left panel) of radius in angular direction and ±1000 km s _1 in the line of sight direction within the volume limited sample (z ^ 0.06 and 
Mr —19.42). In the bottom right panel, the density is measured as a distance to the 5th nearest galaxy in angular direction and ±1000 km s~ 1 in the 
line of sight direction within the volume limited sample. The dotted, short-dashed, solid, long-dashed lines denote all galaxies in the volume limited sample 
(regardless of IRAS detection), VLIRGs, MLIRGs, and LLIRGs, respectively. ULIRGs and HLIRGs are not plotted due to their small number in the volume 
limited sample. 



3.3 Optical Line Properties of Infrared Luminous AGNs 

In this subsection, we investigate optical emission line properties 
of IR luminous AGNs. In Figure 1101 we plot luminosity of the 
[OII1](5008A) emission line against total infrared luminosity for 
AGNs selected in Figure^ It is immediately noticed that there is a 
good correlation between L[oiii] and Lj r , suggesting both L[oiii] 
and L ir can be a good estimator of the total power of the AGNs. 
On the other hand, the excitation line ratios of AGNs do not show 
any correlation with Li r . In the upper panel of Figure llll we show 
[OIII](5008A)/H/9 ratio as a function of L ir . In the lower panel 
of Figure[n] we show [NII](6585A)/Ha ratio against L ir . In both 



panels, there are no good correlation between the excitation ratios 
and Lir. The results suggest that the excitation of AGNs are not 
related to the total power of AGNs. 



3.4 SFR of infrared luminous star-forming galaxies 

In this subsection, we compare SFR estimated from the optical 
emission line (Ha) with Li r . Although it has been often suggested 
that Li r is a good indicator of the galaxy SFR since the infrared 
is thermal radiation from dust heated by the HII regions (Condon 
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Figure 7. Emission line ratios used to select AGNs from our sample. The 
contour shows distribution of all galaxies in the SDSS with r < 17.77 
(regardless of IR detection). The dashed line is the criterion between star- 
bursts and AGNs described in Kewley et al. (2001). Galaxies with line ra- 
tios higher than the dashed line are regarded as AGNs. The magenta, blue, 
green, and red, dots are for ULIRGs, VLIRGs, MLIRGs, and LLIRGs, re- 
spectively. 
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Figure 8. Fractions of AGNs as a function of Li r . The error bars are based 
on Poisson statistics. AGNs are selected using the optical emission line ra- 
tios in FiglTl 
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1992). However, the correlation has never been tested with thou- 
sands of galaxies. 

We measure optical SFR using the prescription given in Hop- 
kins et al. (2003) and Goto (2005a). In order to correct for stel- 
lar absorption, we used US equivalent width (EW) measured in 
Stoughton et al.(2002), and converted it to the stellar absorption 
at the Ha and H/3 wavelengths using the following empirical re- 
lation (Miller & Owen 2002): H/3 EW absorption is equal to US 
EW absorption and Ha EW ais()rl ,t, „ is equal to 1.3 + 0.4 x US 
EW absorption (Keel 1983). Using the stellar absorption corrected 
Ha/H/3 flux ratio, we correct for dust extinction using the prescrip- 
tion given in Hopkins et al. (2003). We apply this extinction cor- 
rection only when both Ha and H/3 lines are detected in emission. 
Finally we apply the 3" fiber aperture correction by scaling the Ha 
flux by the ratio of the Petrosian r-band flux to that measured within 
the 3" fiber of the SDSS (Kewley, Jansen, & Geller 2005). By using 
this stellar absorption, dust extinction, aperture corrected Ha flux 
to the equation given in Kennicutt (1998), we obtain a SFR for each 
galaxy. 

In Fieure fT2l we plot the SFR against Li r . A good correlation 
can be found between the SFR and Li r , suggesting Li r can be used 
as a good indicator of a total SFR of the galaxy. The result also 
supports the idea that in star-forming galaxies, Li r is mainly from 
the re-emission of the light of young stars absorbed by the dust 
(Hopkins et al. 2003; Flores et al. 2004). Note that outliers with 
unusually large SFR are often nearby galaxies with large aperture 
and dust corrections. 

By utilizing the information on Ha/H/3 flux ratio used in com- 
puting the SFR, we can estimate the amount of dust extinction in 
IRGs. In case of CASE B recombination with no dust extinction, 
the theoretical value of the Ha/H/3 ratio is 2.87. We plot Ha/H/3 
ratio against Li r in the lower panel of Figure ^| The solid line 
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Figure 9. Fractions of AGNs as a function of /25//60 IR color. The error 
bars are based on Poisson statistics. AGNs are selected using the optical 
emission line ratios in FiglTl 



connects median values. As can be seen in the figure, there is a 
slight tendency where the Ha/H/3 ratio increase with increasing 
Li r . However, it is important to note that the tendency is small 
compared with the scatter at each Li r , suggesting that there is a 
significant variety in the amount of dust extinction of IRGs. 
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Figure 10. Luminosity of the [OIII](5008A) emission line is plotted against 
total infrared luminosity for AGNs. 



3.5 Optical morphology of infrared luminous galaxies 

In this subsection, we investigate optical morphology of IRGs. As 
is found in FigureQ there are 3438 star-forming IRGs in our sam- 
ple. In addition to the optically measured redshift, all of them have 
CCD-based morphology in 5 optical bands (u,g,r,i, and z) with 
a pixel size of 0.396 arcsec/pixel, providing us with an unprece- 
dented opportunity to investigate optical morphology of statistical 
number of IRGs. 

In Figure 1141 a concentration parameter (Cin) is plotted 
against restframe u—r colour. The concentration parameter is com- 
puted as a ratio of Petrosian 50% flux radius to Petrosian 90% flux 
radius, and is known to correlate well with traditional eye-based 
morphological classification (Shimasaku et al. 2001; Strateva et al. 
2001). It is also known that restframe u — r colour separates early- 
and late-type galaxies well (Strateva et al. 2001 ; Baldry et al. 2004). 
The contours denote all non-AGN galaxies in the SDSS DR3 with 
r < 17.77 regardless of the IRAS detection. As the contour shows, 
early-type galaxies have a peak around (u — r, Cin)={2. 5,0.35). 
Late-type galaxies have a peak around (u — r, Ciri)=(l. 8,0.45). 
The magenta, blue, green, and red dots are for ULIRGs, VLIRGs, 
MLIRGs, and LLIRGs, respectively. AGNs are excluded from this 
plot. 

It is interesting to find many IRGs with Cin ~ 0.35, indi- 
cating significant fractions of IRGs have early-type morphology. In 
the upper panel of Figure [T5l we compare the Cin distribution for 
ULIRGs (dotted line), VLIRGs (short-dashed line), MLIRGs (solid 
line) and LLIRGs (short-dashed line). It is recognized that more 
infrared luminous galaxies have more concentrated or earlier-type 
morphology. This trend is clearer when we plot Cin against Li r in 
the bottom panel of Figure^] The solid line connects median val- 
ues. It is worrisome if difference in average distance between the 
samples artificially created this trend. However, our concentration 
parameter is based on Petrosian parameters, which are independent 
of redshift. And thus, the increase of relative size of PSF to phys- 
ical size of distant object should only make distant objects more 
extended or less concentrated than nearby objects. In Fig. 1151 we 
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Figure 11. Excitation line ratio of AGNs as a function of L; r . The upper 
panel shows [OIII]/H/3 ratio as a function of Li r . The lower panel is for 
[SII]/Ha ratio against Li r . 



found distant objects are more concentrated. This trend cannot be 
caused by the observational bias that more luminous objects are 
at higher redshift. On the contrary, the non-biased relation may be 
even steeper than the observed trend in Fig. 1151 We further compare 
the distribution of Cin in the middle panel of Figure fT5l where we 
plot the distribution of ULIRGs+ VLIRGs with the solid line and 
that of all SDSS galaxies with r < 17.77 with the short-dashed 
line. In Figure fT2l Li r > 10 11 Lq roughly corresponds to the SFR 
of \MQyr~ 1 . For a comparison, we also plot SDSS galaxies with 
the SFR > \MQyr~ 1 with the long-dashed line. Very interestingly, 
although VLIRGs correspond to galaxies with SFR > IMq-ut -1 
(the long-dashed line) in terms of SFR, these two samples have dif- 
ferent distributions of Cin, in the sense that VLIRGs have more 
concentrated, or earlier-type morphology. This difference indicates 
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Figure 12. The SFR computed from the Ha emission line is plotted against 
the total infrared luminosity (£j r ). AGNs are excluded from this figure us- 
ing FigQ 



Figure 14. The concentration parameter is plotted against restframe u — r 
colour. The contours represent all galaxies with r < 17.77 in the SDSS 
DR3. 
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Figure 13. The self-absorption corrected Ha/H/3 ratio is plotted against 
Li r . The solid line connects median values. 



optical spectra are shown in Figure^] Similarly, in Figure [T8l we 
show 20 images of LLIRGs, but LLIRGs are restricted to z > 0.01 
to match the redshift range of images in Figure f!6l The corre- 
sponding optical spectra are shown in Figure fT9l From figs. 131161 
andEU the difference between ULIRGs/VLIRGs and LLIRGs is 
clear: ULIRGs/VLIRGs have frequent merger/interaction or close 
companion, whereas LLIRGs often have normal spiral morphol- 
ogy, suggesting that more IR luminous galaxies are powered by the 
merger/interaction. 

In Figure 1201 we show images of the SDSS galaxies with 
SFR > lM^yr^ 1 and with r < 17.77. These galaxies are se- 
lected to have similar SFR with VLIRGs shown in Fig^| The 
corresponding optical spectra are shown in Figure |2T1 Compared 
with Fig. 1161 although galaxies in both figures have similar SFR of 
> lM0j/r~ , there is a marked difference in their morphological 
appearance: In Fig |20l galaxies are normal spiral or Magellanic- 
cloud-type irregular galaxies, whereas in Fig |16l maioritv of galax- 
ies have signs of merger/interaction. The results cautions in esti- 
mating total SFR of galaxies using Lj r , i.e., if we use Li r to esti- 
mate SFR of galaxies, we will underestimate the SFR of galaxies 
that are not experiencing merger/interaction as is shown in Fig |20l 
Although overall correlation between optically estimated SFR and 
Li r is good (Fig. l 121 . careful attention is needed when Li r is used 
to estimate SFR of individual galaxies. 



that VLIRGs are not just an infrared- view of strongly star-forming 
galaxies, but they might be a different population of galaxies. 

This difference manifests itself more when we look at atlas im- 
ages of IRGs. In Figure|5| we show example g, r, i-composite im- 
ages of ULIRGs. Nine lowest redshift ULIRGs are selected due to 
their larger apparent size on the sky. As is previously reported, ma- 
jority of ULIRGs show dynamically disturbed signatures. The cor- 
responding optical spectra are shown in Figure|4]ln Figure^] we 
show example g, r, i-composite images of VLIRGs. 20 lowest red- 
shift VLIRGs are selected. In addition to ULIRGs, most of VLIRGs 
also show dynamically disturbed morphology. The corresponding 



3.6 Star Formation History of Infrared Luminous Galaxies 

In Figure ^| we notice that there are significant number of 
ULIRGs/VLIRGs with strong B aimer absorption lines such as H<5, 
which motivates us to investigate star formation history of IRGs. In 
Figure l22l we plot restframe u — g colour against H5 EW. The u — g 
colour can be an estimate of the current SFR of galaxies. The H<5 
EW is sensitive to the relative amount of A-type stars in the galaxy, 
and thus, can be an estimate of the past SFR at ~1 Gyr ago. In ad- 
dition, the broadband u — g colour is much easier to be compared 
with the population synthesis models than emission lines which are 
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Figure 16. Examples of g, r, i-composite images of VLIRGs (10 110 Lq < Li r ). The images are sorted from low to high redshift. Only 20 lowest redshift 
galaxies are shown. The corresponding spectra with name and redshift are presented in Fig ll7l (The spectrum of the same galaxy can be found in the same 
column/row panel of Fig |171 ) 



still difficult to be modeled. Therefore, the u — g vs H<5 EW can be 
used to investigate recent star formation history of IRGs. In Figure 
1221 the pentagons, squares, triangles, and circles are for ULIRGs, 
VLIRGs, MLIRGs, and LLIRGs, respectively. The contours denote 
all the SDSS galaxies with r < 17.77. It is recognized that signif- 



icant fraction of ULIRGs/VLIRGs have slightly redder colour in 
u — g than the contours especially at the large H5 EW, indicating 
the star formation history of VLIRGs may be different from that 
of normal galaxies. To clarify this further, we over- plot predictions 
by the population synthesis models (Bruzual & Chariot 2003). The 
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Figure 17. Example spectra of 20 lowest redshift VLIRGs. The spectra are sorted from low redshift. Each spectrum is shifted to the restframe wavelength and 
smoothed using a 20 A box. The corresponding images are shown in Fig. 1161 (The image of the same galaxy can be found in the same column/row panel of 
Fig[l6}. 



solid line is for star forming galaxies, where SFR is constant from 

to 12 Gyr. We use solar metallicity model for simplicity. The dot- 
ted line is for exponentially decaying SFR model with r = lGyr. 
This model tracks the center of the contours, and thus, represents 
a typical star formation history of majority of galaxies. The dashed 
line denotes instantaneous burst model with the burst duration of 

1 Gyr. After the truncation, the galaxy has zero SFR. This model 
is known to reproduce the properties of E+A (or post-starburst) 
galaxies well (Goto et al. 2003c; Goto 2003,2004,2005b). Com- 
pared with these models, MLIRGs and LLIRGs are both consistent 
to be at the end stage of either the exponentially decaying model 
or the truncation model, suggesting that MLIRGs and LLIRGs are 
weakening their SFR in a consistent way as normal galaxies. How- 
ever, ULIRGs/VLIRGs show better agreement with the truncation 
model than with the exponentially decaying model especially at 
high H<5 EW end. This may indicate that ULIRGs/VLIRGs are 
more rapidly reducing their SFR than normal galaxies, at a rate 



close to that of E+A galaxies. Recently Goto (2005) found that 
E+A galaxies have more companion galaxies than normal galax- 
ies, and suggested that the origin of E+A galaxies may be dynam- 
ical merger/interaction with close companion galaxies. In Section 
13.51 we also found that the origin of ULIRGs/VLIRGs may be dy- 
namical merging of multiple galaxies. If the galaxy-galaxy merging 
induces the truncation of star-formation, it is naturally explained 
that both ULIRGs/VLIRGs and E+As are well-reproduced by the 
truncation model in Fig. 1221 However, note that perhaps there is 
no evolutionary link between ULIRGs/VLIRGs and E+A galaxies: 
ULIRGs/VLIRGs have on-going star formation as observed in FIR 
or Ha emission. On the other hand, E+As have no on-going star 
formation by definition. Large H<5 EW of ULIRGs/VLIRGs indi- 
cate they have already started reducing their SFR. Therefore, it is 
not likely that ULIRGs/VLIRGs evolve into E+A galaxies which 
experience a sudden truncation of starburst. 

It is also interesting in terms of evolutionary sequence that we 
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Figure 18. Examples of g, r, i-composite images of LLIRGs (10 10 L Q < L ir ^ 1O 1O ' 5 L ). The images are sorted from low to high redshift. Only 20 
lowest redshift galaxies are shown. The corresponding spectra with name and redshift are presented in Fig. 1191 (The spectrum of the same galaxy can be found 
in the same column/row panel of Fig |19l ) 



obtained a suggestion that ULIRGs/VLIRGs may be reducing their 
star formation activity. This can also be intuitively guessed from 
the presence of strong B aimer absorption lines in ULIRGs/VLIRGs 
CFig. ll7i . In the literature, the evolutionary sequence has been often 
proposed from cool ULIRGs, to warm ULIRGs, to finally QSOs. 



However, in the sample probed in this work, ULIRGs/VLIRGs 
seem to be reducing their star formation activity. If so, perhaps 
ULIRGs/VLIRGs do not evolve into more vigorous starburst galax- 
ies nor powerful AGNs. 
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Figure 19. Example spectra of 20 lowest redshift LLIRGs. The spectra are sorted from low redshift. Each spectrum is shifted to the restframe wavelength and 
smoothed using a 20 A box. The corresponding images are shown in Fig. 1181 (The image of the same galaxy can be found in the same column/row panel of 
FigQU. 



4 DISCUSSION 

In the literature, there has been much evidence that ULIRGs are 
caused by the merger/interaction of multiple galaxies: Allen et al. 

(1985) found 68% of a sample of 19 optically faint IRAS galaxies 
were disturbed, interacting, or had close companions. Soifer et al. 

(1986) found that 12-25% of IRGs with L„ > lO n L are peculiar 
or interacting systems. For 10 ULIRGs, Sanders et al. (1998) con- 
cluded that all were interacting or merging. Lawrence et al. (1989) 
imaged 60 VLIRGs to find 46% are interacting or merging. The 
most luminous ULIRGs were universally classified as strong in- 
teractions/mergers (Sanders et al. 1987,1988b; Hutchings & Neff 
1987: Vader et al. 1987; Kim 1995; Murphy et al. 1996; Clements 
et al. 1996). Sanders et al. (1988) suggested that dissipative galaxy 
merging is important for the formation of ULIRGs. Computation- 
ally, numerical simulations by Bekki (1999,2000) successfully re- 
produced observed Li r of ULIRGs by pair mergers. These mergers 



maybe even multiple mergers of more than two galaxies: Taniguchi 
& Shioya (1998) theorized that ULIRGs might be a merger of more 
than two galaxies (also see Xia & Deng 1997). Bekki et al. (2001) 
showed that such multiple merging of galaxies can trigger repeti- 
tive massive starbursts (SFR ~ lOOMQ-yr^ 1 ) owing to the strong 
tidal disturbance and the resultant gaseous dissipation during the 
merging. Indeed, Borne et al. (1999,2000) observed 99 ULIRGs 
with HST/WFPC2 in /-band to find that 22 of 99 ULIRGs in their 
sample shows evidence for multiple mergers. Colina et al. (2001) 
found that the 67% of the cool ULIRGs observed with HST NIC- 
MOS H have multiple nucleus (also see Cui et al. 2001; Bushouse 
et al. 2002). Our findings of disturbed morphology of ULIRGs in 
Figure [3] are consistent with these results. It is not yet clear how 
much fractions of these merger/interaction is by mergers of a pair 
of galaxies or by multiple galaxies. 

According to our work, this merger/interaction origin of IRGs 
can be applied to lower luminosity VLIRGs. In Section l3~5l we 
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Figure 20. Examples of g, r, i-composite images of star-forming galaxies (SFR > lM@yr~ and r < 17.77). The images are sorted from low to high 
redshift. Only 20 lowest redshift galaxies at z > 0.01 are shown. The corresponding spectra with name and redshift are presented in Fig. 1211 (The spectrum 
of the same galaxy can be found in the same column/row panel of Fig |21l ) 

have found that dynamical merger/interaction between galaxies is 
more frequent in VLIRGs than in LLIRGs (Figs. fT6l and fT5l. It 
is dangerous to search for the signature of merger/interaction by 
eye when one expects to find it. However, in addition to the eye- 
based scrutiny of the morphological appearance of LIRGs, it is 



revealing that this phenomenon also statistically reveals itself in 
terms of the concentration parameter (Cin), which is computed by 
a machine, and thus, is much more objective parameter than the 
eye-based investigation. We found that more IR luminous galaxies 
are more centrally concentrated CFig. 1 15i. This is quite a differ- 
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Figure 21. Example spectra of 20 lowest redshift star-forming galaxies (SFR > IMq^t --1 and r < 17.77) at z > 0.01. The spectra are sorted from low 
redshift. Each spectrum is shifted to the restframe wavelength and smoothed using a 20 A box. The corresponding images are shown in Fig. 1201 (The image of 
the same galaxy can be found in the same column/row panel of Fig |20l . 



ence from less concentrated morphology of normal star-forming 
spiral galaxies, suggesting that IR luminous galaxies is a differ- 
ent population of galaxies from normal star-forming spiral galax- 
ies despite of their similarity in the SFR. In Figs. 1161 and 1201 
there is a distinct difference in morphological appearance between 
VLIRGs and optical selected star-forming galaxies. Although there 
is a good correlation between Li r and optical SFR (Fig |12t . the 
use of L ir to infer the SFR of galaxies is cautioned. Such estima- 
tion may bias high SFR galaxies toward merging/interacting galax- 
ies rather than normal spiral galaxies. The centrally concentrated 
appearance of IRGs is consistent with the hypothesis that IRGs, 
especially ULIRGs/VLIRGs, are created through major merging 
of multiple galaxies: recent major merger computer simulations 
suggest that during collisions, the gas readily loses angular mo- 
mentum due to dynamical friction, decouples from the stars, and 
inflows rapidly toward the merger nuclei (Barnes & Hernquist 
1992,1996; Mihos & Hernquist 1994,1996), creating the central 



starburst, which presumably is observed as high concentration in 
ULIRGs/VLIRGs (Fig.[BJ. In Fig. [22] we found that star forma- 
tion history of ULIRGs/VLIRGs shows better agreement with the 
truncation model. This may also support the merger/interaction sce- 
nario since it has been suggested that in the case of E+A galax- 
ies (Zabludoff et al. 1996; Quintero et al. 2004; Blake et al. 2004; 
Balogh et al. 2004; Goto 2003,2004,2005b) merger/interaction can 
induce sudden truncation of star formation. In Table Q we have 
found that IRGs are not very much brighter than the field galaxies 
in M r . If the origin of VLIRGs is the galaxy-galaxy merger, the 
merging may not be between two very luminous galaxies (3> L*), 
but instead, between less massive sub-L* galaxies (see also Colina 
etal. 2001). 

The environment of IRGs also shows supportive evidence. In 
Section l3~T1 we found that more infrared luminous galaxies exist in 
lower density regions. The result is again consistent with the idea 
that VLIRGs are produced by the merger/interaction of galaxies 
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Figure 15. Normalized distribution of the concentration parameters. In the 
upper panel, we compare the distribution for ULIRG (dotted line), VLIRGs 
(short-dashed line), MLIRGs (solid line), and LLIRGs (short-dashed line). 
In the middle panel, the short-dashed line, the long-dashed line and the solid 
line are for all galaxies with r < 17.77, galaxies with SFR> IMqjjt - , 
and galaxies with Li r > 10 110 Lq. In the bottom panel, Cin is plotted 
against L,; r . The solid line connects the median values. 
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Figure 22. Restframe u — g colour is plotted against H<5 EW. The pen- 
tagons, squares, triangles, and circles are for ULIRGs, VLIRGs, MLIRGs, 
and LLIRGs, respectively. The dashed, solid, and dotted lines are for instan- 
taneous burst, constant star formation, and exponentially decaying models 
(t = lGyr) computed using the GISSEL code (Bruzual & Chariot 2003). 
All models assume solar metallicity for simplicity. The contours denote all 
the SDSS galaxies with r < 17.77. 



since such merger/interaction happens more frequently in the field 
region or in poor groups where relative velocity between galax- 
ies is small. In rich clusters, such merging is less likely to happen 
since relative velocity is too high (~ lOOOfcms^ 1 ; Ostriker 1980; 
Makino & Hut 1997). 

In summary, it has been suggested that the origin of 
ULIRGs/VLIRGs is merger/interaction. In addition to the morpho- 
logical appearance, our work provided statistical evidence such as 
concentration parameter, environment of VLIRGs, and star forma- 
tion history of ULIRGs/VLIRGs to further support the scenario. 



5 CONCLUSIONS 

We have studied optical properties of unprecedentedly large num- 
ber of 4248 infrared galaxies by positionally matching up the SDSS 
and IRAS data. Our findings can be summarized as follows: 

• More infrared luminous galaxies tend to exist in lower density 
regions (Fig.|5J. Especially, VLIRGs have smaller local galaxy den- 
sity than the field galaxies. Since galaxy-galaxy merging is more 
frequent in lower density regions where relative velocity of galax- 
ies is smaller, the trend is consistent with the galaxy-galaxy merger 
origin of VLIRGs. 

• The fractions of AGNs are higher in more IR luminous galax- 
ies (Fig.|8j. 

• There is a good correlation between Lyoui] and Li r , suggest- 
ing both parameters can be a good estimator of the total power of 
AGNs. 

• There is no significant correlation between excitation line ra- 
tios ([OIII](5008A)/H/3 or [NII](6585A)/Ha) and L ir . 

• There is a good correlation between optically estimated SFR 
and Li r , confirming the previous claim that Li r can be a good in- 
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dicator of total SFR for star-forming galaxies. However, there is a 
distinct difference in morphology between optically selected star- 
forming galaxies and VLIRGs, cautioning that the use of IR esti- 
mated SFR may result in a bias toward merging/interacting galaxies 
having higher SFR than normal spiral galaxies. 

• The Ha/H/3 ratio slightly increases with increasing Li r , sug- 
gesting infrared luminous galaxies have more dust extinction. How- 
ever, the scatter is much larger than the trend, suggesting that there 
is a significant variety in the amount of dust extinction of IRGs. 

• More infrared luminous galaxies tend to have more concen- 
trated or earlier-type morphology. This is a contrasting result con- 
sidering that the SFR of VLIRGs corresponds to that of less con- 
centrate spiral galaxies. The result indicates that infrared lumi- 
nous galaxies are not just an infrared- view of strongly star-forming 
galaxies, but there is another factor that causes luminous infrared 
emission. 

• Direct comparison of ULIRGs/VLIRGs and LLIRGs 
images show clear difference: ULIRGs/VLIRGs have more 
merger/interaction or close companion galaxies, whereas majority 
of LLIRGs have normal spiral morphology. It is suggested that the 
merger/interaction with close companions perhaps cause luminous 
infrared emission in ULIRGs/VLIRGs. 

• Comparison with the model SED indicates that 
ULIRGs/VLIRGs may be in a post-starburst phase, being 
consistent with the instantaneous burst/truncation model. Con- 
sidering that the origin of E+A (optically selected post-starburst) 
galaxies is found to be merger/interaction with close companion 
galaxies, this may be additional evidence that ULIRGs/VLIRGs 
may be merging/interacting galaxies. However, ULIRGs/VLIRGs 
are not likely to be the progenitor of E+A galaxies. The timescale 
of truncation of star formation is perhaps much shorter for E+A 
galaxies than ULIRGs/VLIRGs. 
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